Transport Spectroscopy of Single Phosphorus Donors in a Silicon Nanoscale Transistor 
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We have developed nano-scale double-gated field-effect-transistors for the study of electron states 
and transport properties of single deliberately-implanted phosphorus donors. The devices provide 
a high-level of control of key parameters required for potential applications in nanoelectronics. For 
the donors, we resolve transitions corresponding to two charge states successively occupied by spin 
down and spin up electrons. The charging energies and the Lande g-factors are consistent with 
expectations for donors in gated nanostructures. 

PACS numbers: 73.21.-b,61.72.Vv 




The ability to manipulate and measure electrons 
bound to phosphorus (P) donors is a key ingredient 
for the realization of quantum information processing 
schemes using single dopant 8^1. Recent development 
of metal-oxide-semiconductor(MOS)-compatible nanos- 
tructures in silicorPHiS together with single-ion detection 
capabilities^ hold promise for the realization of dopant- 
based spin qubits in silicorP. Such qubits are attractive 
due to the long coherence times of donor electron spins 
in bullP, but further research is necessary to understand 
how the donor spin cohere nce is influenced by the prox- 
imity of a Si-Si02 interfac e 1 10 * 11 ^ . An important step to- 
wards this goal is to create gated nanostructures where 
the energy spectrum of individual donors can be observed 
and, possibly, manipulated by changing the local environ- 
ment or the position of the donors with respect to the Si- 
Si0 2 interface. Pioneering studies of dopant spectroscopy 
have been carried out on a variety of nanostructures22HED 
but the presence and location of the donors were uncon- 
trolled. Here, we describe a double-gated nanoscale field- 
effect-transistor (nanoFET) where a chosen mean num- 
ber of phosphorus donors are deliberately implanted into 
the conduction channel, and use this device to study the 
charge and spin states of these individual P atoms. This 
structure allows convenient control of four key device pa- 
rameters: (i) gate-tunable electron density in the source 
and drain reservoirs^; (ii) number of implanted donors^* 
(iii) depth of implanted donors; and (iv) tunnel coupling 
between donors and reservoirs. In previous studies only 
the tunnel coupling was controlled by choosing the FET 
channel length. 

The nanoFET studied in this paper [see Figure [ljc)] 
consists of two independent aluminum gates: a top gate 
that is biased positively to induce an electron layer that 
constitutes the source and drain reservoirs; and a barrier 
gate that depletes the electron layer in the active region 
of the donors, creating tunnel barriers on each side of the 
donor potential well. This structure enables the electri- 




FIG. 1: (Color online) (a) Coloured SEM image of a tri- 
barrier nanoFET. Positive voltage on the top gate induces 
an electron layer which extends ~25 /xm to source (S) and 
drain (D) ohmic contacts. Barrier gates, which are partially 
underneath the common top gate, allow electrical control of 
the donor energy levels, for example, the chemical potential 
can be brought into resonance with Fermi levels of adjacent 
reservoirs, (b) Zoomed-in SEM image of a single nanoFET 
device. Yellow dashed line indicates the P implant window 
and red dotted line indicates the active region where the mean 
number of donors is three, (c) (Not to scale) Schematic cross 
section of a nanoFET along line XY in panel (b). The red 
source and drain regions are formed by n + diffused contacts. 



2 



cal manipulation of the charge states of the individual 
phosphorus donors under the barrier gate. In particu- 
lar, the barrier gate voltage (Vbg) can be used to tune 
the donor electrochemical potentials into resonance with 
those of the source and drain reservoirs, inducing electron 
transport. The electron density in the source and drain 
reservoirs can be tuned in-situ by the top gate voltage 
(Vtg), while the lithographic width of the barrier gate 
influences the width, and therefore the transparency, of 
the tunnel barriers between donor and reservoirs. The 
excitation spectrum of the donor and its magnetic field 
dependence extracted from the transport measurements 
provide important information on electronic properties of 
donors in close proximity to gate electrodes and induced 
electron layers. 

The devices were fabricated on a high-resistivity 
(~10 kficm) near-intrinsic natural-isotope silicon (100) 
wafei^, with low residual P background doping (^10 12 
cm -3 ). Ohmic contacts for the n + source and drain 
are phosphorus-diffused regions fabricated using stan- 
dard UV-lithography and thermal diffusion at tempera- 
tures ~ 950 °C. The 5 nm gate oxide is high-quality SiC>2 
grown by dry thermal oxidation. A 5 s rapid thermal an- 
neal (RTA) at 1000 °C N 2 ambient was carried out to 
lower the interface trap density to the 2x 10 10 cm _2 eV _1 
range, as measured on similarly processed chipj^. For 
the spatially-selective P implantation, we used a 150 nm 
polymethylmethacrylate (PMMA) layer deposited above 
the SiC>2 as a mask, with 100 nm x 200 nm apertures 
defined using electron beam lithography (EBL). In the 
implantation process, ion acceleration energies of 14 keV 
and 10 keV were utilized and the flux of the P ions was 
controlled to give an average of Nd ~ 3 individual P 
donors in the 50 nm x 30 nm active area of the device, 
that is, the region below the overlap of the top and bar- 
rier gates [red box in FigureQb)]. Subsequently, another 
RTA at 1000 °C for 5 s was applied to activate the donors 
and repair any damage caused by the implantation pro- 
cess. 

The barrier gates used to control the energy levels of 
the P donors and to locally deplete the electron layer 
were fabricated using EBL, thermal evaporation of alu- 
minum, and liftoff. An insulating ~5 nm Al x O y layer 
was then created using plasma oxidation, where the bar- 
rier gates were exposed to a low-pressure oxygen plasma 
(~150 mbar) for 3 min at 150 °C. This Al^Oy provides 
electrical insulation between the barrier gate s and top 
gate for voltage differences greater than ~ 4 V * 5 I 19 I 2 ° I A 
second layer of aluminum forming the top gate was again 
defined using EBL, thermal evaporation, and liftoff. Fig- 
ure [T] shows a device similar to the ones studied here and 
its schematic cross-section. Note that there are three in- 
dependently contacted barrier gates in each device corre- 
sponding to three independent devices, for which individ- 
ual measurements can be done, thus increasing the device 
yield. The top gate is common to all the barrier regions 
and the widening of the top gate between barriers en- 
sures that there are no accidental quantum dots created. 



A useful feature of our structures is that, when prepar- 
ing the apertures for ion implantation, we can choose to 
mask the areas around some of the barrier gates, thus 
obtaining control devices where we know with certainty 
that no P donors are present [see Figure [If a)]. 

We measured 27 P-implanted and 11 control devices 
altogether and present here results from three different 
devices — one non-implanted control device (Sample A), 
and two implanted (Samples B and C) — that represent 
well the qualitative features of the whole batch. Sam- 
ples B and C differ in the energy used for the P ion im- 
plantation (14 keV in B, 10 keV in C) and in the width 
of the barrier gate (30 nm in B, 50 nm in C). All the 
experiments presented here were performed in dilution 
refrigerators at base temperatures T < 100 mK. 

Figure [2jb) shows the source-drain differential conduc- 
tance of Sample B as a function of the dc source-drain 
voltage and the barrier gate voltage, with the top gate 
voltage fixed at Vtg = 3.5 V. We observe the channel 
turn-on as a function of Vbg when the tunnel barrier 
between source and drain becomes transparent. In ad- 
dition, sharp conductance peaks are clearly visible at 
lower barrier voltages below the threshold. These fea- 
tures are signatures of resonant tunneling through dis- 
crete energy levels below the conduction band. Two 
sets of three different peaks at 0.86 V< Vbg <0.93 V 
and 0.96 V< Vbg <l-03 V were observed [Figure [jjfb)]. 
These findings are consistent with the expected num- 
ber of donors (Nd ~ 3) in the active region. However, 
in other devices we often found only one or two sets 
of sub-threshold conductance peaks. This is consistent 
with the Poisson statistics related to the implantation 
processes. In addition, some of the donors can be off- 
center in the active region, and couple too strongly to 
the source (drain) reservoir and too weakly to the drain 
(source) reservoir. The total conductance is then domi- 
nated by the slowest tunnel rate, making certain donors 
immeasurable in transport even if they are physically 
present in the channel. 

We identify peaks a%, b\, and c\ as corresponding to 
the transition from the positively charged state D + to 
the neutral state D° of each donor, whereas peaks b 2 , 
and C2 correspond to a transition to the charged state 
D~ . The peaks associated with the same donor have al- 
most the same lever arm factors a = Cbg/Ce, where 
Cs = C B g + C s + C D + C T g + Co and C B g> C T g, Cg, 
and Cd are the coupling capacitances between the cor- 
responding donor and the barrier gate, top gate, source, 
and drain, respectively. The self-capacitance of the donor 
is denoted by Cq. This a-factor can be determined from 
the positive and negative slopes of each conductance peak 
in the differential conductance plot [Figure [2Tb)] , where 
the positive slope is given by Cbg/(Ce — Cs) and the 
negative slope equals Cbg/Cs (with drain grounded)^. 

Figure 2(b) shows that the slopes of peaks b\ and b 2 are 
similar to each other, with a — 0.45, and hence they are 
attributed to the same donor. Multiplying a by the dif- 
ference in barrier gate voltage A Vbg = 62 mV between b\ 
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FIG. 2: (Color online) (a) Stability diagram showing differen- 
tial conductance as a function of barrier gate and dc source- 
drain (SD) bias of a control device (Sample A) with no im- 
planted donors, (b) Stability diagram of Sample B with six 
labeled peaks, (c) A cut of the stability diagram in panel 
(b) at Vsd = V, highlighting the resonant tunneling peaks 
due to three donors, (d) Direction of peak shifts under global 
magnetic field up to 7 T. Peak a\ (red) shifted with the field 
to the left while peak &2 (green) shifted to the right. Schemat- 
ics indicate that ai is due to tunneling of spin-down electrons 
and 62 is a two-electron singlet state. All data in this figure 
were measured at T < 100 mK , and Go = e 2 /h. 



and 62 yields the D° to D~ charging energy for donor b, 
E c = 28±3 meW^I Peaks 0,2 and C2 are just barely resolv- 
able in Figure [2jc) , so we cannot make a specific match 
with cither of the D° peaks a\ and c\. However, we can 
use the slope a = 0.26 of peak 01, together with the 
gate voltage difference A Vbg = 131 mV between peaks 
ai and (02,02), to calculate the charging energy for donor 
a to be E c = 34±6 meV. In both cases E c is lower than 



the 44 meV expected for phosphorus in bulk ! 23 * 24 ] due to 
the fact that there is a substantial capa citive coupling of 
the donor to the surrounding electrode d 12 * 13 l Across the 
entire series of implanted devices we have observed E c 
ranging from 21 meV to 39 meV, with an average value 
of 28 meV. A simple model of the gate electrodes and 
the donor at the mean implantation depth in silicorP^U 
provides a rough estimate of the capacitances required 
to calculate the charging energy E c = e 2 /Cs. With a 
donor depth d = 5-20 nm below the Si-SiO*2 interface, 
the model yields E c = 39-40 meV. The discrepancy with 
the measured values can be, for example, due to the limi- 
tations of a simple capacitance model for this ultra-small 
few-electron system. 

Magnetospectroscopic studies of Sample B are shown 
in Figure [2jci) , obtained by applying a magnetic field up 
to 7 T in the plane of the electron layer. As the mag- 
netic field was increased, peaks ai and c\ shifted to lower 
values of Vbg whereas the peak 6 2 shifted to higher val- 
ues^. The direction of these shifts provides information 
about the spin-polarity, relative to the external field, of 
the electrons that contribute to the tunneling current: a 
shift to lower (higher) Vbg corresponds to tunneling of 
spin-down (-up) electrons. The magnetic-field-dependent 
shifts are consistent with the assignment of peaks a\ and 
ci to D° and 62 to D~ states, as inferred from the bias 
spectroscopy in Figure [2^b) . This is due to the fact that 
electrons in D~ , being a two-electron ground state, form 
a singlet with opposing electron spins^. A quantita- 
tive analysis of the magnetic field dependence is given 
for Sample C below. 

As a control, another nanoFET structure (Sample A) 
was fabricated on the same substrate but with no im- 
planted dopants in the active region. The stability dia- 
gram is shown in Figure [2]ja), with Vtg = 2.0 V. It is 
significantly different from the data obtained for devices 
with donors in the active region. In particular, there are 
no sharp resonant tunneling features with large charging 
energies well below threshold. A weak minimum is ob- 
served in the threshold region (1.32 V< Vbg <l-40 V) 
which is most likely due to the formation of an unin- 
tentional dot, due to weak local disorder in the barrier 
region. This is consistent with the small charging energy 
(<10 meV) and the large width of the peaks, which im- 
plies transparent tunnel barriers. Features reminiscent 
of Kondo-enhanced conduct ance^HUl i n this region, plus 
the absence of other peaks at lower Vbg, suggest that 
this could be the first electron in such a dot. The con- 
trol data reinforce the interpretation that the peaks ob- 
served in Sample B are due to resonant tunneling through 
the deliberately implanted P donors. We note that some 
control devices processed without a final forming gas an- 
neal (FGA) showed sub-threshold conductance features 
which could be attributed to localized states associated 
with interface traps, for example. These features were 
clearly distinguishable from donors due to their much 
lower charging energies, typically only up to 13 meV. We 
found that a FGA can reduce significantly the charge 
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trap density as well as the charge noised 
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FIG. 3: (Color online) (a) Bias spectroscopy showing dif- 
ferential conductance as a function of barrier gate and dc 
SD-bias of Sample C. (b)-(e) Detailed bias spectroscopies of 
peak d 1 at B = 0, 2, 4, and 6 T. Red (yellow) arrows indicate 
spin-down (-up) states entering the bias window, (f) Energy 
splitting (AEz) of the spin- up and spin-down states deduced 
from panels (b)-(e). The dashed line indicates 2/isB/e for 
reference, where /j,b is the Bohr magneton, (g) Positions of 
the peaks di (red circles) and c?2 (green circles) in the bar- 
rier gate voltage as a function of global magnetic field. The 
dashed lines are linear fits to the data points. 

Sample C was fabricated using the same processes 
as Sample B, with the exception that the phosphorus 
donors were implanted with 10 keV ion-beam accelera- 
tion energy, and the barrier gate was 50 nm wide, yield- 
ing weaker tunnel coupling between each donor and the 
reservoirs, and accordingly sharper conductance features. 
With 10 keV implantation energy the donors are ex- 



pected to be ~12 nm from the Si-SiC>2 interface^, thus 
closer to the gate compared with the donors implanted 
with 14 keV beam energy in Sample B. The bias spec- 
troscopy of Sample C is presented in Figure §[a). Here, 
only one pair of D and D~ states, denoted by d\ and d 2 , 
is strongly visible although the mean number of donors 
under the active region (here 50 x 50 nm 2 ) is N4 ~ 5. 
This is not surprising because the wider tunnel barriers, 
as compared with Sample B, make the requirement of the 
donors being well-centered even more stringent!^. 

In the Coulomb blockaded region between d± and d 2 , 
additional faint conductance features were observed at 
finite source-drain bias. These features are attributed to 
nearby weakly coupled donors or other charge traps^. 
We also observed random telegraph signals (RTS) in the 
conductance in well-defined regions at | Vsd I > 10 mV. 
This phenomenon most likely arises because the electro- 
chemical potential (and thus the conductance) of the vis- 
ible donor depends on the charge state of the surrounding 
donors or traps, due to mutual charging energy effects. 
Random switching of weakly coupled charge centers can 
then give rise to the observed RTS. Thus the donor 
can be possibly used as a charge sensor for the nearby 
charge centers^. Our preliminary theoretical analysis 
shows that the distance between the visible donor and 
the nearby donor or charge trap is more than 25 nmP. 

The sharpness of the conductance peaks in the wide- 
barrier Sample C allowed us to study in more detail the 
spin states of the donor. Figures 3(b-e) show a set of 
bias spectroscopy data on the peak di of Sample C, with 
applied magnetic fields B = 0, 2, 4, 6 T. For B > 2 T we 
clearly resolve an excited state at finite bias, which we 
attribute to the excited Zeeman state (spin-up) of the 
first electron on the donor. The peaks in the differen- 
tial conductance due to spin-down and -up electrons are 
indicated by red and yellow arrows, respectively. The 
separation of these peaks on the Vbg axis multiplied by 
a should then equal the Zeeman splitting of the D° spin 
states. To verify this, we fitted Gaussian functions to the 
differential conductance peaks in Figures 3(b-e) at vari- 
ous fixed SD bias points to obtain an accurate measure 
of their position. The resulting energy splitting AEz 
is plotted as a function of magnetic field in Figure |3jf ) , 
together with error bars obtained from the fitting proce- 
dure and the error of a (including AEz at B = 0). As 
shown in Figure |3][f), AEz is in good agreement with the 
Zeeman splitting for a Lande g-factor g = 2, as expected 
for an electron bound to a hydrogenic donor potential. 

Figure [3^g) shows the Zeeman shifts of peaks d\ and 
o?2 in barrier gate voltage. As expected, the peaks shifted 
linearly in opposite directions, but the magnitudes of the 
shifts are different. The peak d\ shifted an average of 
-24 /iVT -1 whereas the peak d 2 shifted an average of 95 
/xVT -1 , yielding individual g- factors far from 2. How- 
ever, if the difference of the shifts for d\ and d 2 is con- 
sidered, we obtain the correct value g — 2.4 ± 0.5. The 
asymmetry in the peak shifts can be caused, for exam- 
ple, by a slow global drift in the electric potentials in the 
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device over time. 

While the nanoFET devices in this study were used 
to investigate phosphorus donors, we note that these 
double-gated structures could be used to study discrete 
energy states of any donor or acceptor atom in silicon. 
In the latter case, this could be achieved using a pMOS 
nanoFET, employing p + source-drain contacts and a 
negative gate voltage. By employing on-chip waveguides 
or microwave cavities to create a magnetic held in res- 
onance with the spin split D° donor state, one could 
use these structures for the study of electrically-detected 
magnetic resonance of single donor spins, as recently 
demonstrated for small ensembles^. 

In conclusion, we have shown that resonant tunneling 
of electrons through discrete energy states of phospho- 
rus donors in silicon can be observed in planar nanoscale 
FET structures, where the presence and depth of the 
donors is deliberately chosen, and basic device param- 
eters such as tunnel couplings and electron density can 
be conveniently controlled. We have identified the D° 
and D~ charge states of individual donors, measured 
their charging energy, and observed the Zeeman shift of 
their spin states in magnetic helds up to 7 T. Control de- 
vices, deliberately fabricated with no donors in the chan- 
nel, showed markedly different transport behavior than 



the characteristic features present in donor-implanted de- 
vices. These experiments provide important understand- 
ing of the electric and magnetic properties of individual 
phosphorus donors needed for the realization of Si:P na- 
noelectronics. Moreover, these str uctu res are compatible 
with single-ion implant detector d 7 ^ 37 ^ thus enabling de- 
vices with precisely one implanted donor. The demon- 
strated ability to resolve spin states of donors in these 
devices could also enable the measurement of electrically- 
detected magnetic resonance for a single-donor electron 
spin. 
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